LPTENS 09- 



Critical interface: twisting spin glasses at Tc 



E. Brezin") , S. Franz^) , G.Parisi 



Laboratoire de Physique Theorique, Ecole Normale Superieure 
24 rue Lhomond 75231, Paris Cedex 05, France, e-mail: brezin@lpt.ens. 
Laboratoire de Physique Theorique et Modeles Statistiques 
Universite Paris-Sud 11, Centre scientifique d'Orsay, 15 rue Georges 
Clemenceau 91405 Orsay cedex, France 
'^^ Dipartimento di Fisica, "Sapienza" Universita di Roma, P.le A. Moro 2, 

00185 Roma, Italy 2 
INFM-CNR SMC, INFN, "Sapienza" Universita di Roma, P.le A. Moro 2, 

00185 Roma, Italy 



Abstract We consider identical copies of spin glasses in finite dimension 
coupled at the boundaries. This allows to identify the spin glass analogous 
of twisted boundary conditions in ferromagnetic system and leads to the 
definition of an interface free-energy that is positively defined and that should 
scale with a positive power of the system size in the spin glass phase. In 
this note we study the behavior of the interface at the spin glass critical 
temperature within mean field theory. We show that the leading scaling 
of the interface free-energy does not depend on replica symmetry breaking, 
and can be obtained by simple scaling arguments using a cubic theory for 
critical spin glasses. 
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1 Introduction 



Sensitivity to boundary conditions is a fundamental tool to study the nature 
of the Gibbs states of extended physical systems. If ergodicity is broken 
and Gibbs states are not unique, different boundary conditions can select 
pure phases or induce interfaces in the system. For example, in Ising-like 
ferromagnets in the ferromagnetic phase, the choice of homogeneous "plus" 
(respectively "minus") boundary conditions where all the spins outside a 
large region are fixed to point in the up (resp. down) direction is enough to 
select the pure state with positive (resp. negative) magnetization. Twisted 
boundary conditions, where along a specific direction + and — conditions are 
chosen at the opposite boundaries (while neutral boundary condition, e.g. 
periodic ones, are chosen in the other directions) induce an interface with 
positive tension. On the contrary, in the paramagnetic phase the system 
is insensitive to boundaries and the surface tension is equal to zero. Using 
boundary conditions one can therefore investigate the stability of low temper- 
ature phases and the lower critical dimension below which the ferromagnetic 
phase is unstable. Around the critical temperature, the behavior of the inter- 
face tension is described by finite size scaling, which implies that right at the 
critical temperature, the interface tension vanishes, and the free-energy cost 
to impose twisted boundary conditions remains finite in the thermodynamic 
limit. 

In trying to extend this kind of considerations to study the stability of 
spin glass phases, one is confronted with the fact pure phases are in this case 
glassy random states strongly correlated with the quenched disorder. They 
cannot therefore be selected using boundary conditions uncorrelated with 
the quenched disorder. To project on some pure phases boundary conditions 
should be chosen that depend self-consistently on the quenched disorder. 
One way to do that, is to consider two "clones" of a spin glass with identical 
disorder coupled at the boundaries. In ref. [10] it was considered a situa- 
tion where different values of the overlap (a measure of correlation between 
configurations) was imposed on the boundaries along one direction. In that 
paper it was investigate the stability of low temperature phase replica sym- 
metry broken (RSB) phases below the critical temperature, against spatial 
fluctuations of the overlap that restore replica symmetry, and it was found 
a value of the lower critical dimension below which the fluctuations destroy 
RSB equal to D^cd = 5/2. 

In the present note we we consider systems exactly at the critical tem- 
perature, where the techniques developed in [lOj do not apply. 

In the case of ferromagnetic systems the behavior of the interface free- 
energy at the critical temperature is a well understood issue in critical fluctu- 
ation theory [T]. One has that below dimension 4 the interfacial free-energy 
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scales with the system size as cr ~ /,-(<^-i) while it scales as a ~ above 
dimension 4. 

We then compare two kind of boundary conditions. We consider two 
identical copies of cylindrical spin glasses (i.e. we impose p.b.c. in the trans- 
verse directions) and we couple them on the boundaries. In a first kind of 
b.c. the systems are chosen to be identical on both boundaries. In a second 
kind of b.c, the two systems are still imposed to be identical on one of the 
boundaries, but they are opposite in the other boundary. 

We have already remarked that in ref. [lOj it was considered a situation 
where different values of the overlap where fixed at the boundary. In this 
situation the interface energy is a a function of the temperature, of the ge- 
ometrical parameters and of the values of the overlap. In ref. [TT] it was 
showed that it may be interesting to consider the extreme case where the 
overlap is 1 at one boundary and -1 (or 1) at the other boundary. These ex- 
treme choice simplifies the analysis and it can be implemented in numerical 
simulations in a very simple way. 

The plan of the paper is the following: In section II we discuss in detail 
the Ising ordered case. This is well understood and it will guide us in the 
analysis of the more complex spin glass case. Section III which constitute 
the core of this paper is devoted to the analysis of the Edwards- Anderson 
model at T^. Finally we draw our conclusions. 



2 Pure system 

We consider an Ising-like system in a cylindrical geometry. The d-dimensional 
cylinder has a length L and a cross-section area A ~ L*^"^. The boundary 
conditions are periodic in the directions transverse to the axis of the cylinder, 
but on the surfaces at the end of the longitudinal direction the spins may be 
all up or all down. Below Tc there is an interfacial tension a defined as the 
difference of free energy per unit area 

a = (1) 

It is a function of the temperature t ~ (T — Tc) , and of the aspect ratio 
L/A^^^~^. This tension vanishes near Tc as 

ar-^i-tr. (2) 

A scaling law due to Widom [1] relates fi to the correlation length exponent 

^ = iy{d-l) (3) 
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and it may be derived in the standard renormahzation group framework [2] 
(for d<A). 

Assume that we want to know the behavior of this tension for finite L at 
Tc. We use finite size scahng 

ait,L) = {-trf{L/0 (4) 

and, since a is non singular at Tc for L finite, it imphes that 

fix) ~.^o x-^/\ (5) 

Therefore at Tc the interfacial tension vanishes with 1/L as 

a(0,L) ~L-^/" = L-^"^-^). (6) 

In other words at Tc 

Tt,^ - T^,t ~ L'. (7) 

In appendix A we give the details of the verification of this behavior ([7]) 
in a simple Landau theory (i.e. tree level of a (f'^ theory) : i.e. for a free 
energy at Tc 



(we have assumed translation invariance in the {d — 1) transverse directions). 
Plus-plus b.c. are defined be the fact that at both boundaries the order 
parameter is fixed to the value m, plus- minus conditions correspond to impose 
the value m on the left boundary and — m on the right one. One finds that 
independently of the values of m, the surface tension behaves as 

with a positive constant C given by 

C=^-^rS(-) ^r6(i)~44.8. (10) 

1927r2 127r3/2 W ^ ' 

Clearly this mean field theory should only be compared with ([7]) in four 
dimensions, the result (Q is indeed in agreement since A ~ T^. 



3 Spin glass 

Twisting fixed boundary conditions uncorrelated with the disordered cou- 
plings for a spin glass does not produce any interface since it can be gauged 
away in a change of the random couplings. It may increase or decrease the 
free energy, depending on the random instance.. 

In this paper, we would like to consider two copies of a similar cylindrical 
sample of spin glass with identical couplings Jij and we assume that we fix 
two successive boundary conditions as follows : 
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• In the up-up situation the overlap between the spins of the two copies 
in the left {x = 0) plane and in the x = L plane are equal to a positive 
value Iq. 

• In the up-down situation the overlap is still equal Iq in the left plane, 
but it is —Iq in the right plane. 

In ref. [IT] Contucci et al. considered boundary conditions of this kind 
in the extreme case Iq = 1. Notice that in this case, the construction leads 
to effective interactions within the spins in the planes x = and x = L 
are doubled. Choosing the interactions in these planes with the same statis- 
tics as the others would lead to stronger effective interactions between the 
spins in these planes then all the other interactions. In order to avoid such 
strong difference, in [11], the strength of the interaction among spins on the 
boundary planes was reduced by half. 

We aim here at characterizing the behavior of the difference AF = 
■^T.-i^-^T.T g^|- fuuctlou of L. Agalu we shall limit ourselves to mean field 

theory. We submit to test the Widom scaling law, verifying that above the 
upper critical dimension Djjcd = 6 the interface tension scales as cr ~ L^^. 

We would like to study this problem within a Landau-like small order 
parameter expansion of the free-energy close to T^. Unfortunately, this can 
not be done directly if the value Iq of the order parameter at the boundary 
is chosen to be large. In particular, we could not directly consider the con- 
ditions chosen in [iTj. However, the small order parameter expansion can 
be saved even in this case, arguing that, analogously to the ferromagnetic 
case, the result for AF should not depend on the value |/o| imposed at the 
boundaries. We will see that this is true within the Landau model. We have 
also confirmed this in a spherical spin-glass model [5], which has a Landau 
expansion identical to the one of the standard Ising spin glass, and for which 
one can write equations valid even if for large values of the order parameter. 

We start from a short range standard Edwards- Anderson model. A After 
tracing out over the spins we can rewrite it in the long wave length limit in 
terms of a 2n x 2n matrix Q (for the bulk we need only a n x n matrix Q) 
when a and b are both smaller than n, Qab refers to the overlap between the 
spins in the first copy ; when the two indexes are larger than n we are dealing 
with the overlaps Qa,/3 in the second copy and for one index lower than n 
and one larger , Qaa is a cross-overlap. The critical mean field free-energy is 
here 

F = J d'xlhr (VQ)^ + f tr Q3 + I ^ Q^,] (11) 

ab 

Assuming again translation invar iance along the transverse directions Q{x) = 
Q(z), we assume that at the replica-symmetry is unbroken. Therefore we 
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assume 

Qab = Qa(3 = q{z) (12) 

and 

Qaa = k{z){l - 6a,a) + l{z)6a,a- (13) 

Let us note that the free energy ffTTl) is invariant under S'„ x S'„ , i.e. inde- 
pendent permutations of the n rephcas of system one, and of the n rephcas 
of system two. However the last term in ( |T3l) breaks this symmetry down to 
Sn- Indeed one may imagine the boundary conditions as follows : one adds 
a coupling Jar between the spins of the two systems located in the plane 
z = L. The boundary conditions a = t would correspond to J going to +00, 
whereas a = — r is enforced by the limit J — — 00. After replicating the 
coupling becomes J J2 <^oJa and thus the boundary conditions lead to a ba,a 
term in (IT51) [6j 
Then one has 

tr (^)2 = 2n(n - 1)^'^ + in(n - 1)A;'2 + InV" (14) 
dz 

tr = 2n{n - \)(n - 2)q^ + Qn{n - l)(n - 2)qk'^ + I2n{n - l)qkl (15) 
Ql, = n{n - l)q' + ^n{n - l)k' + nl' (16) 
and thus, in the zero replica limit, the free energy reads 

^ = Af dz[-q'^ - k'^ + 1'^ + |(V + 12gfc2 - I2qkl) + ^ (-g^ -\^^ + ^ 

(17) 

The quartic term may be dropped at Tc since the functions q{z),k{z),l{z) 
remain small for large L. The equations of motion are thus simply 



2q" + Aw{q'^ + P - kl) =0 
2k" + Aw{2qk - ql) =0 

2/" + Awqk = (18) 

The plus-plus, and plus-minus boundary conditions are different in the 
two cases for the functions l{z) and k{z). In the plus-plus case, if one imposes 
the same boundary conditions on q and k , q{z) = k{z) is solution and 
we are left with two equations. There is one constant of motion, namely 
4g'^ — /'^ -|- y(8g^ — 3lq), but the last integration has to be done numerically. 
In the plus-minus case, one has to deal with three different functions. 

We impose the boundary conditions in the following way: we consider a 
sample of size (1 + 2r)L with —rL < x < (1 + r)L and fix the value of l{x) to 
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preassigned values on the regions —rL < x < 0, that we call left boundary, 
and L < a; < (1 + r)L, that we call right boundary. In the case of plus-plus, 
the value Zq > is imposed both right and left boundaries, in the plus-minus 
case the value Zq > is chosen on the left boundary and —Iq is chosen on the 
right boundary. The value of q{x) and k{x) on the boundaries are not fixed 
by the constraint and should correspond to extremization of the free-energy. 
Since r > and L is large, g(0) and fc(0) and q{L) and k{L) should take the 
same values that that they would take if a uniform overlap profile l{x) = ±/o 
for all X was chosen, namely g(0) = k{0) = q{L) = k{L) =1^/2 for plus-plus 
conditions and g(0) = k{Q) = q{L) = lo/2, k{L) = —lo/2 for plus-minus 
conditions. 

We expect that, as verified explicitly in the case of the ferromagnet, the 
interface free-energy does not depend on the imposed value Iq. In that case, 
a scaling argument based on the fact that the dominant interaction terms 
in the free-energy are the cubic ones (namely the transformation x —>■ xL 
q —>■ w~^L~'^q and analogously for A; an /), suggests that the interface tension, 
in absence of accidental cancellations should behave as AL~^ for large L, 
i.e. for A ~ L^~^, AF ~ L^~^. Indeed, as it is easy to see by the above 
rescaling, we can write F++ = ++{10101"^), F+_ = -^g__{lowL'^). This 
is of course true only if no accidental cancellations appear. If this is the case, 
we are indeed authorized to neglect the quartic term. Since this is the term 
responsible for RSB, we can neglect RSB effects to the leading order. 

In order to verify that — remains finite for large L, we have 

integrated numerically the rescaled equations f|T8|) . for various values of the 
parameter v = IqwL"^. This is done by simple a relaxation method, where 
initial profiles are iterated until one observes convergence to a solution of a 
discretized version the equations ( |T8|) . In fig. ([21) we plot the function f{v) = 
g+-(v) — g+^{v), together with a power law fit of the kind f{v) = a + 
While f{v) continue to have an appreciable dependence on v (and thus on L) 
for very large values of v, the fit indicates that foo = lim^,_,oo /(f) is different 
from zero. We have in addition explicitly verified that the inclusion of the 
quartic term in the free-energy gives a subleading contribution to AF. 

4 Conclusions 

In this paper we have considered the interface free-energy to impose different 
boundary conditions in a spin glass. We have concentrated in the behavior 
at the critical point and found that the behavior of the interface follows 
the usual pattern of critical phenomena. The free-energy difference above 6 
dimensions is of order L^~^ as naive scaling suggests. This implies that RSB 
effects can be neglected to the leading order at the critical temperature. This 
is very different from what found at low temperature [TU] , where the scaling 
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Figure 1: The order parameters in the two solutions. Upper paneh the 
functions l{x)/lo and q{x)/lo = k{x)/lo as a function of x/L for v = 1000 
and r = 1/2. Lower Paneh the functions l{x)/lo and q{x)/lo and k{x)/lQ for 
the same values of the parameters. 



of the interface was found to depend critically on zero modes associated to 
RSB. 
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Figure 2: The rescaled free-energy difference, / = {F^ — F^^)w'^L^^^ as a 

function of f = IqwL'^, together with a power law fit of the form f{v) = a+^. 
Chi-square fitting for v > 10^ gives the parameters a = 480±20, b = 1806±15 
and c = 0.19 ±0.01. 



Appendix A : Pure systems 

The free energy is 



r+L/2 _ I 
-L/2 

1. Up- up boundary conditions 

The magnetization m in the planes z = ±L/2 is given and we solve the 
equation of motion 

- 9?" +uip^ = (20) 

with the b.c. 

= ipi+L/2) = m (21) 

The mechanical analog is a negative energy bounce off the potential 
— j(p^ starting ai (p = m bouncing at v?o and returning to ip = m. It is 
given by 

Mz) dtp fu,^ , , , , 

^"^ ^^P^ - V 2 

Then (pQ is determined by 

r^^ = -f-L (23) 



For L large the order parameter (/9o at the centre of the sample is small, 
much smaller than the given finite value m on the boundaries, and it 
is thus given asymptotically by 



-j^'^^ + 0{l/L) (24) 



'8 dt 
The free energy is then given by 

Jm 4 

The corrections to the first leading term are of order We finally 

note that 

r * =^M) (26) 

and 

.oo rft ^ r^(i/4) _ 1 .271 

2. Up-down boundary conditions 

We arc still dealing with the mechanical analog of a motion in the 
inverted potential —Uip'^/A but, now with a positive energy solution 
going from ip{L/2) = m to ip{—L/2) = —m. The solution is given by 

Jm V2^2 ^ ^ ' 

and (pi is determined by 

r^^^f-L (29) 

i-m ^^4 + ^4 V 2 

Again (^i is of order 1/L and given asymptotically by 



Then the free energy 

^ta ^ _|^^4 ^ (31) 
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The last asymptotic estimate 



dtVt^=-{—f+ (32) 



(33) 
(34) 



/o 3Vi Jo y^WTT + f^ 

and the values of the integrals 

^ dt _ r2(l/4) 

dt T\l/A) 



lo VW+T + t^ V36^ 
completes the calculation. 

We may now compute F-^^i — Fj^-f ; the constant term cancels and we 
are left with a difference proportional to A/L^ as expected 

^T,i-^T,T_ Ir 5 ^8^^^^^ L^r6(l/4)] (35) 



Appendix B : A spherical Spin Glass Model 

As we observed in the main text, the boundary conditions defined in ref. 
imply high overlaps on the boundaries, and cannot be analyzed in the 
context of Landau small order parameter expansions of the free-energy. We 
define here a model with the following properties: 

• The replica free-energy of the model admits a simple closed form in 
terms of the overlap matrix Qa,b{x), and it can be explicitly analytically 
continued to n — > when the Parisi ansatz is assumed. 

• The Landau expansion close to Tc coincides up to the fourth order term 
with the one of the SK model. 

Consider a system where on each site, labelled by its longitudinal and trans- 
verse coordinates x and y respectively, there are K spins cr^ y, r = 1, ...,K 
subjected to the spherical constraint Z]r(>S'^,y)^ = K. We write the Hamil- 
tonian of the model as a sum of terms that couple spins on the same plane 
and terms that couple spins in adjacent planes. 

H = j2nr + j:nr (36) 

x=0 x=0 

where the Hamiltonians Ti"^^ and Ti^"'^ are gaussian random variables with 
variances 

iKMnrir]) = E/^(^ E & 

y ^ zeVy 



(37) 
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where we denoted by Sy the value of the spin on site y of the x-th plane, 
^y!z ~ 'W ^r=i ^x,y'^w,zi oveilaps between spin configurations on different 
sites. The different functions of the overlap are chose to be /(g) = ^{q'^+yq^), 
fx{q) = f{q) for a: 7^ 0, L and /(O) = \f{q) for x = 0, L. Notice that with 
this choice it is possible to express Hamiltonian (1361) in terms of two body 
and four body gaussian couplings between the spins. We can now introduce 
two copies of the system with ++ and H — boundary conditions. The ++ 
conditions consist in considering two copies o"^ y and y constrained to be 
identical for x = 0, L: o"q y = Tq y and (j£ y = r]^ y. The +— conditions consist 
in considering two copies y and r^y constrained to be identical for x = 
but opposed for x = L: = r^y and cr£y = — t£ y. As for the reduced 
model, the replica treatment of the problem involves the introduction of 
two local n X n overlap matrices, which under the replica symmetric ansatz 
and the assumption of independence of the overlap profiles of the transverse 
spatial coordinate, can be parametrized in terms of the functions q{x), l{x) 
and k{x) x = 0, ...,L of the main text. The resulting F++ free-energies and 
as a function of these parameters can be decomposed as 

/, k}] = F'"'"'[{q, I, k}] + Fq""'^^ + F^^"'^"'^?^ (38) 
F+_ [{q, I, k}] = ^'"''^[{g, /, k}] + Fo""'^^ + Ft°_""'^"^^ (39) 

where 

_f^pbuik ^ p2j2[f{l) + /(/(a;)) - f{q{x)) - f{k{x))] 

x=l 

+(3' J2 /(I) + /((^(^) + + l))/2) - fM^) + Qi^ + l))/2) - fm^) + + l))/2) 

x=l 

1 7/ N , \ 7 / M g(x) + A;(x) 

+- Log 1 + / X - g X - X + p . 
^ 1 + l[x) - g(x) - k[x) 

+Log [1 - /(x) - g(x) + k{x)] + '/^f (40) 

1 — l{x) — g(x) + k[x) 



(41) 



1 

+ 2 



Log [1 - g(0)] + + Log [1 - g(L)] + ^^^^ 



l-g(O) l-g(L) 



P pboundary 



++ 

32/ 



+ /((I + 1{L - l))/2) - /((g(L - 1) + g(L))/2) - /((g(L) + k{L - l))/2) 

(42) 
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boundary 

+ /((I + KL - l))/2) - /((g(L - 1) + g(L))/2) - /((g(L) - k{L - l))/2) 

(43) 

Given the complexity of the expression, we have manipulated them through 
the Mathematica software in order to obtain the equations of motion, and 
integrated the resulting equations numerically by the relaxation method. The 
integration gives result fully compatible with the analysis of the reduced 
model, confirming the independence on the detailed boundary conditions 
imposed. As a proxy of the free-energy difference — in figure [3] we 
plot the difference in free-energy density in the center of the box multiplied 
by L^, which according to the argument given in the main text should tend 
to a constant for large L. As in the case of the reduced model, fitting curves 
of the kind jiV) = a + h/ L'^ confirm this behavior. 




0.1 I ' • — • — • • — • — • ' ' — I 

1 10 100 1000 

L 

Figure 3: The difference of free-energy density in the center of the box for the 
two kinds of boundary conditions as a function of the system size L. Notice 
that even and odd values of L give rise to different curves. The curves are 
plotted together with a power law fit of the form f{L) = a + j^. Chi-square 
fitting for L > 20 gives the parameters a = 0.15 ± 0.04, b = 3.1 ± 0.3 and 
c = 0.44±0.04 for L even and a = 1.3±0.1, b = 46.9±1.7 and c = 0.58±0.04 
for L odd. 
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